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Die f olgenden Angaben aind den vom Anmelder eingereichten Unterlagen entnommen 

® Signalauswerteverfahren zur Detektion von QRS-Komplexen in Elektrokardiogramm-Signalen 

@ Ein Signalauswerteverfahren zur Detektion von QRS- 
Komplexen in Elektrokardiogramm-Signalen weist fol- 
gende Verfahrensschritte auf: 

- Abtasten des EKG-Signals (4) und Umwandlung in dis- 
krete, zeitlich aufeinanderfolgende Signalwerte (x(n), 
Xfg(n)), 

- Vergleichen der Signalwerte (x(n), Xf q (n)) mit einer adap- 
tiv aus diesen ermittelten Schwellenwertfunktion (K(n)), 

- Feststellen einer Haufigkeitszahl (D(n}) innerhalb eines 
definierten Segmentes der aufeinanderfolgenden Signal- 
werte (x(n), Xf q (n)), mit der die Betrage der Signalwerte 
(x(n), Xfq(n)) die Schwellenwertfunktion (K(n)) unter- 
schreiten, und 

- Vergleichen der ermittelten Haufigkeitszahl (D(n)) mit ei- 
nem definierten Grenzwert (GW), wobei ein Unterschrei- 
ten des Grenzwerts (GW) signifikant fur das Vorliegen ei- 
nes QRS-Komplexes (5, 6, 7) in dem definierten Segment 
des EKG-Signals (4) ist. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Signalauswerteverfahren zur Detektion von QRS-Komplexen in Elektrokardio- 
gramm-(EKG- und IEGM-)Signalen. 

5 (0002] Zum Hintcrgrund der Erfindung ist festzuhalten, daB die automatische Analyse von EKG-Signalen zur Perfek- 
tionierung der Funktionsfahigkeit von Herzschrittmachern und Defibrillatoren eine immer groBere Rolle spielt. Solche 
kardiologischen Gerateimplantate j lingerer Bauart bieten dementsprechend auch die Fahigkeit einer EKG-Analyse. In 
diesem Zusammenhang hat die Erkennung von QRS-Komplexen und R-Zacken in EKG-Signalen einen uberragenden 
Stellenwert. Diese Bedeutung resultiert aus der vielfaltigen Verwendbarkeit der Information iiber das zeitliche Auftreten 

10 des QRS-Komplexes, z. B. bei der Untersuchung der Heraaten variability, bei der Klassifikation und bei der Datenkorn- 
prcssion sowic als Basissignal fur sckundarc Anwcndungcn. Ubcrhaupt nicht odcr falsch dctckticrtc QRS-Komplcxc und 
R-Zacken sind dabei problematisch im HinbLick auf die Leistungsfahigkeit der auf die Detektion folgenden Verarbei- 
tungs- und Analysestufen. 

[0003] Einen breiten Oberblick iiber bekannte Signalauswerteverfahren zur Detektion von QRS-Komplexen in EKG- 
15 Signalen bietet der Fachaufsatz von Friesen et al. "A Comparison of the Noise Sensitivity on Nine QRS Detection Algo- 
rithms" in IEEE Transaction on Biomedical Engineering, Vol. 37, Nr. 1, Januar 1990, Seiten 85-98. Die dort vorgestell- 
ten Signalauswerte-Algorithmen basieren durchgehend auf einer Auswertung der Amplitude, der ersten Ableitung des 
Signals sowie dessen zweiter Ableitung. Bei den vorgestellten Algorithmen wird unterschieden zwischen solchen, die 
eine Analyse der Amplitude und der ersten Ableitung vornehmen, solchen, die nur die erste Ableitung bewerten und sol- 
20 chen, die die erste und zweite Ableitung in Betracht Ziehen. Kurz umrissen wird bei alien Algorithmen abgepruft, ob der 
jeweilige Signalparameter vorgegebene Schwellen uber- oder unterschreitet, wonach bei Auftreten eines solchen Ereig- 
nisses mit einem vorgegebenen Schema das Auftreten weiterer definierter Ereignisse abgepruft und bei ErfUUung be- 
stimuiter Kriterien auf das Vorliegen eines QRS-Komplexes geschlossen wird. 

[0004] Ein anderer Aspekt bei der Signalauswertung zur Detektion von QRS-Komplexen ist bei der Implementierung 
25 solcher Verfahren in implantierten Herzschrittmachern zu beachten. Im Hinblick auf die naturgemaB bei solchen Geraten 
gegebenen Beschrankungen bei der Energieversorgung und der Rechenkapazitat ist darauf zu achten, daB zur Detektion 
von QRS-Komplexen moglichst einfache Algorithmen mit moglichst wenigen Rechenoperationen auf der Basis von 
ganzen Zahlen anstatt von reellen Zahlen durchzufuhren sind. 

[0005] Ferner wurden zur QRS-Erkennung Signalverarbeitungsmethoden aus den Gebieten der linearen und nicht- li- 
30 nearen Filterung, der Wavelet-Transformation, der kiinstlichen neuronalen Netze und der genetischen Algorithmen an- 
gewendet. Bei hohenSignal-Stor-Abstanden und nicht pathologischen Signalen, also im Falle guter Signaiverhaltnisse, 
erzielten diese A uswerte verfahren zuverlassige Ergebnisse. Waren keine solchen Verhaltnisse gegeben, konnte die Lei- 
stungsfahigkeit der Aus werte verfahren unter Umstanden drastisch absinken, was im Hinblick auf eine zuverlassige 
Funktionsweise eines Schrittmachers naturlich nicht akzeptabel ist. 
35 [0006] SchlieBlich ist aus der alteren, jedoch nachveroffentlichten deutschen Patentanmeldung 100 11 733.3 der An- 
melderin eine QRS-Erkennung auf der Basis von Nulldurchgangszahlungen bekannt. Dabei werden folgende Verfah- 
rensschritte durchgefuhrt: 

- Abtasten des Signals und Urnwandlung in diskrete, zeitlich aufcinanderfolgende Signalwerte, 
40 - HochpaGfillerung der abgetastelen Signalwerte, 

- Bestimmung des Vorzeichens jedes Signalwerts, 

- laufendes Abpriifen der Vorzeichen aufeinanderfolgender Signalwerte auf Vorliegen eines Nulldurchgangs zwi- 
schen zwei aufeinanderfolgenden Signalwerten, 

- Ermittlung der Anzahl der Nulldurchgangc in cincm definicrten Segment der aufeinanderfolgenden Signalwerte, 
45 und 

Vergleichen der ermittelten Nulldurchgangszahl mit einem definierten Schwellen wert, wobei ein Unterschreiten 
des SchweUenwertes signifikant fur das Vorliegen eines QRS-Komplexes in dem definierten Segment des Signal- 
verlaufes ist. 

50 [0007] Urn dabei im Bereich auBerhalb der QRS-Komplexe eine signifikant hohe Nulldurchgangszahl zu erhalten, 
wird dort zu dem Hoch- oder BandpaBgefilterten und quadrierten EKG-Signal ein hochfrequentes Uberlagerungssignal 
mit im Vergleich zur Amplitude des QRS-Komplexes niedriger Amplitude addiert. 

[0008] Diese Vorgehensweise steht. der eingangs erlauterten Forderung nach moglichst. einfachen Algorithmen natur- 
lich entgegen. 

55 [0009] Ausgehend von den geschilderten Problemauken liegt der Erfindung die Aufgabe zugrunde, ein Signalauswer- 
teverfahren zur Detekdon von QRS-Komplexen in EKG-Signalen anzugeben, das mit vergleichsweise geringer Rechen- 
kapazitat und auch bei problemalischen Signalverhaltnissen einsetzbar ist dabei jedoch zuverlassige Delekuonsergeb- 
nisse liefert. 

[0010] Diese Aufgabe wird durch die im Patentanspruch 1 angegebenen kennzeichnenden Verfalirensschritte wie folgt 
60 gelost: 

- Abtasten des EKG-Signals und Urnwandlung in diskrete, zeitlich aufeinanderfolgende Signalwerte, 
Vergleichen der Signalwerte mit einer adaptiv aus diesen ermittelten Schwellenwertfuiktion, 

- Feststellen einer Haufigkeitszahl innerhalb eines definierten Segmcntes der aufeinanderfolgenden Signalwerte, 
65 mit der vorzugsweise die Betrage der Signalwerte die Schwellenwertfunktion unterschreiten, und 

- Vergleichen der ermittelten Haufigkeitszahl mit einem definierten Grenzwert, wobei ein Unterschreiten des 
Grenzwerts signifikant fur das Vorliegen eines QRS-Komplexes in dem definierten Segment des EKG-Signals ist. 
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[0011] Kernpunkt des erfindungsgemaBen Verfahrens ist die Anwendung eines Schwellenwertvergleichs und einer 
nachfolgenden Haufigkeitszahlung, wobei die Morphologie des QRS-Komplexes ausgenutzt wird. Der QRS-Komplex 
irn EKG-Signal ist namlich durch eine relaliv hochamplitudige Oszillation charakterisiert, die den Signal verlauf ausge- 
pragt von der in der Regel vcrrauschten und mil einem Offset beaufschlagten Basis-Linie des Elektrokardiogramms weg- 
fUhrt. Die Frequenz dieser kurzen Oszillation lieet in einem Bereich, in dem andere Signalkomponenten, wie die P- und s 
T- Welle, nur geringen EinfluB ausiiben und vorzugsweise durch Vorfilterung - z. B. eine Hoch- oder BandpaBfilterung - 
entfernt werden konnen. Nach der Unterdriickung dieser tieffrequenten SignaLkomponenten ergeben sich aufgrund von 
hdherfrequentem Rauschen Signalschwankungen urn die Basis-Linie, die in dem Bereich, wo kein QRS-Komplex auf- 
tritt, das EKG-Signal dominieren. Der QRS-Komplex erscheint dann in diesem Signalkontext als langsame, nur kurz an- 
dauernde Schwingung hoher Amplitude, die von der Basis-Linie des EKG-Signals signifikant wegfuhrt. Werden also die 10 
Signalwcrtc mit cincr das Signalrauschcn rcprascnticrcndcn Schwcllcnwcrtftunktion vcrglichcn, so sind die Betrage der 
Signalwerte auBerhalb des QRS-Komplexes meist kleiner als diese Schwellenwertfunktion. Insoweit ist die Haufigkeits- 
zahl, mit der die Betrage der Signalwerte die Schwellenwertfunktion unterschreiten, groB. Im Bereich des QRS-Komple- 
xes iiberschreiten nun die Betrage der Signalwerte die Schwellenwertfunktion, und zwar signifikant. Es stellt sich also 
eine geringe Haufigkeitszahl fur Unterschreitungen der Schwelle im Zuge des QRS-Komplexes ein. Durch ein Vergiei- 15 
chen der ermittelten Haufigkeitszahl mit einem definierten Grenzwert kann also der QRS-Komplex selektiert werden. 
Ein Unterschreiten des Grenzwertes signalisiert namlich das fur den QRS-Komplex typische Uberschreiten der Schwel- 
lenwertfunktion. 

[0012] Das erfindungsgemaBe Verfahren der QRS-Komplex-Erkennung hat. sich als robust, gegenuher Rauschstorun- 
gen und rechentechnisch besonders einfach realisierbar erwiesen. Insoweit eignet es sich besonders fur die Implementie- 20 
rung bei der Realzeit- Analyse von EKG-Signal-Morphologien in Herzschrittmachem. 

[0013] Die bereits erwahnte HochpaBfilterung wird vorzugsweise mit einer unteren DurchlaBfrequenz von 18 Hz vor- 
genoimnen. Damil konnen die tieflrequenten Komponenten, wie die P- und T- Welle sowie eine Grundliniendrift unler- 
driickt werden. Ferner wird der QRS-Komplex dadurch zur Signalkomponente mit der niedrigsten Frequenz, die das Si- 
gnal wahrend ihres Auftretens dominiert 25 
[0014] Zur VergroBerung des Signal-Stor-Abstandes kann femer vorgesehen sein, die Signalwerte vor dem Verglei- 
chen mit der Schwellenwertfunlktion und der Haufigkeitszahl zu quadrieren. Damit werden kleine Signalbetrage relativ 
zu groBeren Signalbetragen abgeschwacht, was die Detektierbarkeit des QRS-Komplexes weiter verbessert. 
[0015] Der Wert der Schwellenwertfunktion wird vorzugsweise adaptiv aus einer flieBenden Mittelung der BandpaB- 
gefilterten und quadrierten Signalbetrage bestimmt. 30 
[0016] Im folgenden wird das erfindungsgemaBe Verfahren anhand der beigefugten Zeichnungen in einem Ausfuh- 
rungsbeispiel naher erlautert. Es zeigen: 

[0017] Fig. 1 eine hochst schematische Darstellung des Signalverlaufs eines QRS-Komplexes in einem EKG-Signal, 
und 

[0018] Fig. 2 ein Strukturschema des erfindungsgemaBen Signalauswerteverfahrens zur Detektion von QRS-Komple- 35 
xen in EKG-Signalen. 

[0019] Wie aus Fig. 1 deutlich wird, besteht ein schematisch dargestellter QRS-Komplex aus einer relativ hochampli- 
tudigen Oszillation urn die Null-Linie 1, die das EKG-Signal 4 in der Q-Zacke 6 zuerst in negativer Richtung liber die 
Null-Linie 1 wegfahrt. AnschlieBend wird das EKG-Signal 4 in der R-Zacke 5 mit einem steilen Anstieg in den positiven 
Bereich und einem anschlieBenden steilen Abfall wiederum in den negativen Bereich unter Ausbildung der S-Zacke 7 40 
ubergefuhrt. 

[0020] Real ist das EKG-Signal 4 mit einem bestimmten Pegel verrauscht, was in Fig. 1 durch den strichlierten Signal- 
verlauf angedeutet ist. Wird nun dieses verrauschte Signal 4 abgetastet und in diskrete, zeitlich aufeinanderfolgende Si- 
gnalwcrtc umgcwandclt und BandpaB-gcfiltcru so konnen dicsc Signalwcrtc mit der in Fig. 1 schematisch als mit Krcuz- 
chen untcrlcgtc Linic dargcstclltcn Schwellenwertfunktion K(n) vcrglichcn werden. Wie sich nun aus Fig. 1 anschaulich 45 
und modellhaft ableiten laBt, liegt im Bereich auBerhalb des QRS-Komplexes der Betrag des EKG-Signals meist unter- 
halb dieser Schwellenwertfunktion K(n). Es ergeben sich also beispielsweise im Bereich Nl signifikant hohe Haufigkei- 
ten fur Signalbetrage lx(n)l unterhalb der Schwellenwertfunktion K(n). 

[0021] Im Bereich der Q-Zacke 6 weicht nun der Betrag des EKG-Signals sehr stark von der Schwellenwertfunktion 
K(n) in positive Richtung ab. Die Haufigkeitszahl D(n) innerhalb des Segmentes N2 des QRS-Komplexes 6 riir dieses 50 
Ereignis ist also wesenttich geringer als die Haufigkeitszahl innerhalb des Segments Nl. Insoweit kann die Haufigkeits- 
zahl D(n) fur die Erkennung des QRS-Komplexes herangezogen werden. Dieser wird dann als vorliegend erkannt, wenn 
die Haufigkeitsanzahl D(n) einen definierten Grenzwert 0 unterschreitet. 

[0022] Als ein Kernpunkt der Erfindung gegenuber dem Stand der Technik ist an dieser Stelle hervorzuheben, daB auf- 
grund der Feststellung der erorterten Haufigkeitszahl und deren Vergleich mit einem definierten Grenzwert nicht - wie 55 
im Stand der Technik - die Amplitude des EKG-Signals daraufhin uberprUft wird, ob ein besdmmter Schwellenwert ab- 
solul uberschritlen wird, um daraus auf den QRS-Komplex zu schlieBen. Viehuehr wird gewissennaBen uberpruft, wie 
lange das EKG-Signal eindeutig auf einer fur das Vorliegen eines QRS-Komplexes sprechenden Seite einer Schwellen- 
wertfunkdon "verharrt". Erst das Vorliegen einer bestimmten Zeitdauer dieses Zustandes wird zum RuckschluB auf das 
Vorliegen eines QRS-Komplexes herangezogen. Demzufolge werden nur kurze, starke MeBschwankungen nicht als (fal- 60 
sche) QRS-Komplexe (sogenannte "False positive"-Fehler) erkannt. 

[0023] Der detaillierte Ablauf des erfindungsgemaBen Auswerteverfahrens wird nun im folgenden anhand der Fig. 2 
naher erlautert. 

[0024] Das EKG-Signal 4 wird abgetastet und in diskrete, zeitlich aufeinanderfolgende Signalwerte x(n) umgewandelt. 
Die Abtastrate betragt beispielsweise f T = 360 Hz, d. h„ daB das EKG-Signal in eine Folge von 360 MeBwerten pro Se- 65 
kunde umgewandelt wird. Das abgetastete EKG-Signal x(n) wird dann eingangsseitig einer BandpaBfilterung BP unter- 
zogen, die der Entiernung aller Signalkomponenten, die nicht zum QRS-Komplex gehoren, dient. Dazu zahlen die P- 
und T- Welle sowie hochfrequcntere Storungcn, die z. B. von der bioelektrischen Muskclaktivitat ausgehen. Der einge- 



3 



DE 101 05 431 A 1 

setzte Filter BP ist linearphasig, nicht rekursiv unci besitzt eine BandpaBcharakteristik mil den DurchlaBfrequenzen f gl = 
1 8 Hz und f g 2 = 27 Hz sowie den Sperrgrenzfrequenzen f sl = 2 Hz und f s o = 50 Hz. Die Filterordnung ist FO = 26. Die 
Gruppenlaufceitdes BandpaBfilters BPentspricht dementsprechend 13 Abtastwerten und muB bei der Bestimmung des 
Zeitpunktes des QRS-Kornplexes beriicksichtigt werden. 
5 [0025] Die so erhaltenen Signalwerte xKn) werden anschlieBend in einer Quadrierstufe QS gem^B folgender Bezie- 
hung quadriert: 

IxKn)! 2 



10 [0026] Die so aus den urspriinglichen Signalwerten x(n) durch eirie Art Betragsbildung aufbereiteten Werte Xf q (n) wer- 
den nun cincm Vcrglcichcrkomplcx HZ zugcfuhrt, dcr dicsc Signalwerte mit cincr adaptiv aus dicscn crmittcltcn Schwcl- 
lenwertfunktion K(n) vergleicht. In Fig. 2 ist der die Ermittlung der Schwellenwertfunktion K(n) betreffende Verfahrens- 
komplex mit AS gekennzeichnet. Darin wird ein geeigneter Wert fur die Funktionskoeffizienten K(n) adaptiv aus den Si- 
gnalwerten Xfq(n) geschatzt. Dazu werden flieBend die BandpaB-gefilterten und quadrierten Signalwerte mit Hilfe eines 

15 Gedachtnisfaktors 

(0 < A* < 1) rekursiv gemittelt. 
K(n) = A* K (n - 1) + ( 1 - X k ) x fq (n) ■ c 

wobei c ein Konstante ist. 
20 [0027] Empirisch ergeben sich als geeignete Werte A* = 0,98 und c = 8. 

[0028] Die durch den GedachtnisfaktorA* vorgegebene Mittelungszeit besiimmt im wesentiichen die Adaptionsge- 
schwindigkeit dieser Schatzung, wobei sowohl zu kurze als auch zu lange Mittelungssegmente die Leistungsfahigkeit 
des Signalaus werte verfahrens beeinlrachtigen konnen. 

[0029] Im Verfahrenskomplex HZ werden nun die Signalwerte Xf q (n) mit der Schwellenwertfunktion K(n) - wie ange- 
25 sprochen - verglichen. Dabei wird festgestellt, in welcher Richtung die Signalwerte Xf q (n) von der Schwellenwerifunk- 
tion K(n) abweichen. Daraus wird eine Ilaufigkeitsanzahl D(n) innerhalb dieses definierten Segmentes N festgestellt, die 
die Anzahl oder Haufigkeit von Ereignissen reprasenuert, bei denen die Betrage der Signalwerte Xf q (n) unterhalb der 
Schwellenwertfunktion liegen. D(n) kann auch auf rechengunsuge Weise rekursiv uber 
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D(n) = X D D(n - I) + (I - X D ) d(n) mitd(n) =< 



0furx fq (n) 2k(n) 
1 furx fq (n) <k(n) 



ermittelt werden. Ein kleiner Wert fur die Haufigkeitszahl D(n) spricht dafur, dafi der Beu*ag des EKG-Signals 4 dauer- 
haft oberhalb der Schwellenwertfunktion K(n) liegt, was einen verlaBlichen Parameter fur das Vorliegen des QRS-Korn- 
plexes darstellt. 

[0030] Im Zuge des eriindungsgemaBen Verfahrens muB also noch ein Grenzwert 0 ermittelt werden, dessen Unter- 
40 schreiten signifikant fur das Vorliegen eines QRS-Komplexes in dem definierten Segment des EKG-Signals 4 spricht. 
0(n) wird rekursiv aus D(n) berechnet nut 



0(n) = X* 0(n - 1) + (1 - hd D(n) 

45 wobei ein Gcdachtnisfaktor 0 < X$ < 1 verwendet wird. Fur das am Endc dieser Bcschrcibung kurz wicdcrgcgcbcnc Zah- 
lenbeispiel wurde dieser Gedachtnisfaktor zu \$ = 0,99 gewahlt. Liegt D(n) unterhalb des Grenzwertes 0, ist ein QRS- 
Komplex erfaBt worden, andemfalls nicht. 

[0031] Die Abpriifung, ob die vorstehende Bedingung erfiillt worden ist, findet in der Entscheiderstufe Es gemaB Fig. 2 
statt. 

50 

Patentanspruche 



I. Signal auswerteverf ah ren zur Detektion von QRS-Komplexen in Elektrokardiogramm-(EKG-)Signalen, ge- 
kennzeichnet durch folgende Verfahrensschritte: 
55 - Abtasten des EKG-Signals (4) und Umwandlung in diskrete, zeitlich aufeinanderfolgende Signalwerte 

<x(n)), 

- Vergleichen der Signalwerte (xf<n), Xf q (n)) nut einer adaptiv aus diesen ermiltellen Schwellenwertfunktion 
(K(n)), 

- Feststellen einer Haufigkeitszahl (D(n)) innerhalb eines definierten Segmentes der aufeinanderfolgenden Si- 
60 gnalwerte (Xf<n), Xf q (n)), mit der die Signalwerte (xKn), Xf q (n)) die Schwellenwertfunktion (K(n)) unterschrei- 

ten, und 

- Vergleichen der ermittelten Haufigkeitszahl (D(n)) mit einem definierten Grenzwert (0), wobei ein Unter- 
schreiten des Grenzwertes (0) significant fur das Vorliegen eines QRS-Komplexes (5, 6, 7) in dem definierten 
Segment des EKG-Signals (4) ist. 

65 2. S ignalaus werte verf ah ren nach Anspruch 1 , dadurch gekennzeichnet, daB die EKG-Signalwerte (x(n)) nach dem 

Abtasten einer HochpaBfilterung, vorzugsweise einer BandpaBfilterung (BP) unterworfen werden. 
3. Signal auswerteverf ah ren nach Anspruch 2, dadurch gekennzeichnet, daB die untere und obere DurchlaBgrenz- 
frequenz (f gl , f g2 ) des BandpaBfilters (BP) bei ca. 18 Hz und ca. 27 Hz liegen. 
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4. Signalauswerteverfahren nach einem der Anspruche 1 bis 3, dadurch gekennzeichnet, daB die Signalwerte (x(n)) 
einer Betragsbildung unterworfen werden. 

5. Signalauswerteverfahren nach Ansprach 3, dadurch gekennzeichnet, daB die vSignalwerie (Xf<n)) vor dern Ver- 
gleichen mit der Schwelienwertfunktion (K(n)) und Ermitteln der Haufigkeilszahl (D(n)) zur Betragsbildung qua- 
driert werden. 

6. Signalauswerteverfahren nach einem der Ansprilche 1 bis 4, dadurch gekennzeichnet, daB der Wert der Schwel- 
ienwertfunktion (K(n)) adaptiv aus einer flieBenden Mittelung der Signalwerte (x(n)), vorzugsweise der BandpaB- 
gefilterten und quadrierten Signalwerte (xf q (n)) ; uber eine durch den Gedachtnisfaktor (Xk) vorgegebene Mitte- 
lungsperiode (P) bestimmt wird. 

7. Signalauswerteverfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, daB der fur einen QRS- 
Komplcx (5, 6, 7) signifikantc Grcnzwcrt (0) der Haufigkcitszahl (D(n)) als adaptive Schwcllc aus der Haufigkcits- 
zahl (D(n)) selbst variabel eingesteilt wird. 
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